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ABSTRACT 


A method to measure the size of particles In LDV applications Is dis- 
cussed. Since In LDV the velocity of the flow Is associated with the velocity 
of particles to establish how well they follow the flow, in the present 
method the interferometric probe volume is surrounded by a larger beam of 
different polarization or wavelength. The particle size is then measured from 
the absolute intensity scattered from the large beam by particles crossing the 
fringes. Experiments using polystyrene particles between 1.1 and 3.3 pm and 
larger glass beads are reported. It is shown that the method has an excellent 
size resolution and its accuracy is better than 10% for the particle size studied. 


NOMENCLATURE 

b Q waist radius 

d particle diameter 

G gain function of instrument 

I intensity 

K. scattering cross-section 

n particle relative index of refraction 

r' parallel polarization 

S perpendicular 

x,y,z coordinates 

V particle visibility 


Greek Symbols 
a = Trd/X size parameter 

y intersection angle of small laser beams 

0 angle between laser beam and axis of collecting lens 

p laser wavelength 

0 solid angle of collection 
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Subscripts 


1 refers to small beams 

2 refers to large beam 

3 center of beam 

INTRODUCTION 


A technique to measure the size of particles in laser Doppler velocimetry 
(LDV) systems is described here. In typical LDV applications the velocity of 
the flow is inferred from the measurement of the velocity of very small 
particles seeded in the fluid. These particles are chosen to be small (0.5 pm 
to 1 u®) and when possible neutrally bouyant to insure that their velocity is 
the same as the velocity of the flow. There are, however, many situations of 
interest where large particles might be generated, and since they don't 
necessarily follow the flow they will introduce an error in the velocity 
measurement. Examples of these situations are flows with very large accelera- 
tions, sooting environments, wind tunnels with large amounts of oil droplets, 
and in general, two phase flows where it is of interest to distinguish the 
velocity of the two phases. Many LDV electronic processors have built-in 
saturation circuits to reject signals with excessive amplitudes. This, 
however, will not prevent large particles from being measured. Since the 
intensity profile of the laser beams is typically Gaussian, large particles 
moving through the edge of the beams will scatter light with an acceptable 
signal level. Therefore, the saturation circuits mentioned above will only 
accept large particles when they cross through the edge of the Gaussian laser 
beans. This will not eliminate velocity bias but it might introduce addi- 
tional errors since the large particles are detected over a different volume 
which could have a different velocity. The method discussed here eliminates 
the above limitation and provides for a mechanism to obtain the size distribu- 
tion of the measured particles. In its simplest form this method can be used 
co produce a velocity histogram of small particles and one of large particles 
where the separation can be arbitrarily selected. 

In this method two small laser beams of a given wavelength and polariza- 
tion are crossed in the middle of a large laser beam of the same wavelength 
but normal polarization. The two small beams interfere where they cross, and 
this interference pattern identifies the center of the large beam. This 
defines a region of almost uniform intensity within the large beam where the 
light scattered by the individual particles can be related to their size. 

Only particles crossing through the interference pattern are measured and 
their velocity is obtained from the Doppler signal. Their size is obtained 
from the absolute scattered light since the incident intensity is known. 

The method which is referred to as Polarization Intensity Maximum (PIMAX) 
was described in Reference (1) where it was used to measure the size and 
velocity of droplets in a spray. For LDV applications It is necessary to show 
that the PIMAX technique can also discriminate the size of particles down to a 



fraction of a micron. Numerical and experimental studies were conducted to 
demonstrate such ability and they are reported here. 

The method could axso be Implemented In two-color LDV systems. Here two 
small beams of one color cross in the middle of two crossing large beams of 
another color forming two fringe patterns perpendicular to each other. An 
before, the small beams identify the certter of the large beams from which the 
size information Is obtained. 

The model presented here Is for spherical particles with a diameter of 
the same order as the wavelength of the incident beams. Therefore, the full 
Mie solution (2) is necessary to analyze the scattering patterns. Results are 
presented for flows seeded with known size particles. 


DESCRIPTION OF OPTICAL TECHNIQUE 

The method discussed here bases the size measurement on the absolute 
intensity scattered by the particle crossing the probe volume, and the 
velocity measurement on the classical Doppler signal. A similar method could 
be described for a single beam veloclmeter but the approach used here is for 
LDV. In situ single particle counters are limited because of the nonuniform 
profile (typically Gaussian) of laser beams. Under this condition a particle 
crossing the middle of the beam will scatter more light than a similar 
particle crossing through the edge. Therefore, the relationship between size 
and scattered light is not unique. 

To circumvent this problem, two small beams of a given wavelength or 
polarization are crossed in the middle of a larger beam of different wave- 
length or polarization identifying a region of almost uniform intensity and, 
therefore, removing the Gaussian ambiguity. The crossing beams will interfere 
and a fringe pattern will be formed in the middle of the large beam 
(Figure 1). Signals exhibiting a sinusoidal modulation correspond to parti- 
cles that cross the fringe pattern and therefore the middle of the large 
beam. Both size and velocity of individual particles can be extracted from 
this signal. The method and results reported here correspond to the two- 
polarization (PIMAX) system. 

If we refer to the small beam as 1 and the large be »s 2, the intensity 
profiles in the probe volume can be separated by their t • ization and given 
by: 


21 e *P (-j^-)[x 2 +y 2 +z 2 Y 2 /4] 
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where I Q Is the center intensity, y Is the Intersection angle, b fl the waist 
radius, X the laser wavelength, and x,y,z the coordinates. The z dependence 
of the large beam Is negligible. If we also assume that ^ » 0 (which is an 
excellent assumption since a pinhole In the receiver will limit the value of 
z), the Intensity scattered by a spherical particle Is given by: 


I s - 2 I q K 1 (d,n,0,Q,X,,P)G 1 exp[(- -2_)(x Z +y 2 )][l + cos 2 If* . v] (3) 
11 b 

°1 



- 2_)( x 2 +y2j] (4) 
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where K is the scattering cross-secticn as obtained by the Lorentz-Mie theory 
(3). It Is a function of a complex size function d, the index of refraction 
n, the collection angle 8, the solid angle of collection a, the wavelength X 
and the polarization. Here it is assumed that the small beams have a polari- 
zation parallel to the scattering plane while the polarization of the big beam 
i3 perpendicular. G is the gain function of the instrument, and V is the 
visibility of the measured particle. 

It should be pointed out that the visibility of the particle is not an 
adequate parameter to obtain the size, given that many of the particles of 
interest are below 3 pm, and in general, many of the particles will be sub- 
micron. 

Equation (3) gives the signal response of the laser Doppler veloclmeter 
which will establish the detectability of the signal. The processing logic 
will be the following: signals xhlbltlng sinusoidal modulation will have 
crossed the fringe pattern which is located in the middle of the large beam. 
Therefore, signals validated by the laser veloclmeter correspond to particles 
crossing the middle of the large beam and their scattered light can be 
inverted to size since the Gaussian ambiguity is removed. In Equation (4) the 
exponent can therefore be approximated as unity and we get 


1 - I K 2 (d,n,0,fl,X,S)G 2 exp [( 




(5) 


The K 2 coefficients will be obtained from Lorentz-Mie scattering using a modi- 
fied version of the numerical program developed by Dave (4). Equation (5) can 
then be solved for the particle size which is contained in K£ (d,n,6,ft,X,S). 
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NUMERICAL COMPUTATIONS OF THE SCATTERING FUNCTIONS 


To obtain the particle diameter, d, It is necessary to know the 
functional relationship of K 2 (d,n, 6,0, X,S) as given by Equation (A). This 
function can be quite complex and It is necessary to find the conditions under 
which the ambiguities, if any, are within tolerable error margins. The compu- 
tations were made in a PDP11 computer. Parametric studies were conducted to 
establish optimum experimental conditions. These parameters include the angle 
of collection ( 0 ), the solid angle of collection ( 0 ) , and index of refraction 
(nj - Im( 02 )). The results show that the best optical characteristics are 
near forward scattering angles of collection (0 < 5°). 

Since the index of refraction of the particles present in IDV applica- 
tions may be quite different, it is important to establish conditions which 
are less sensitive to these variations. Both real (n^) and imaginary ( 02 ) 
parts of the refractive indices were varied to check the sensitivity of these 
parameters. Figures 2 and 3 show the scattered intensity as a function of the 
particle size parameter a for different values of n^ and n 2 at different 
scattered angles. 

It can be concluded from these calculations that shallow angles of 
collection (< 5°) offer the most favorable conditions. 

The calibration of the system was conducted with known size polystyrene 
particles, therefore the scattering function corresponding to these particles 
was also evaluated. The results are shown on Figure A and, as before, the 
shallow angles offer the best conditions. 


APPARATUS AND EXPERIMENTAL FACILITY 

Figure 5 shows a schematic of the optical system used In the experiments. 
The system consists of a transmitter and a receiver positioned 5° off axis. 

The transmitter uses a 7 mW He-Ne laser. The laser beam is focused on a 
diffraction grating (DG) by lens Lj, and the three major orders are collimated 
by lens L 2 . The zero order beam goes through a beam compressor formed by L^ 
and L 5 and its polarization is rotated by PR. Lens L 3 focuses and crosses all 
three beams to form a probe volume like the one shown in Figure 1. The light 
scattered by particles crossing the probe volume is collected by L^ and 
focused ento the photomultipliers by L 7 . A cube polarizer (CP) efficiently 
divideb the scattered light into its two polarization components. In 
addition, a polarization filter was placed in front of the FMT which looks at 
the big beam (BPMT) to reduce the crosstalk between the two polarizations. 

The size of the prohe volume was typically 600 pm for the large beam and 
86 pm for the small ones. The fringe spacing was 10.7 pm. To reduce the size 
of the probe volume we used a slit in front of the BPMT. The width of this 
slit was typically equal to the diameter of the small probe volume (100 pm) 
and its length was about four times longer to avoid signal masking. Thus, the 
probe volume was about 10“5 cm" 3 which in general is quite adequate for LDV 
applications . 
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The outputs of the photomultipliers are Input to an electronic breadboard 
which measures the Doppler period and the peak Intensity of each scattering 
center. A description of the electronics can be found in Reference (5). 


THE PARTICLE GENERATOR 

The particle size of interest is that corresponding to typical LDV 
applications. We used polystyrene particles of 1.1 pm, 1.7 pm, 2.7 pm and 
3.3 pm in diameter. These are latex particles made by Dow Chemicals of good 
size uniformity and of spherical shape. The particles come suspended in water 
with a concentration of 10% by weight. A few droplets of the particle suspen- 
sion were introduced and diluted in a nebulizer. An air compressor provided 
the air flow to produce a mist carrying the polystyrene particles out of the 
nebulizer and into a heated chamber. There the water was evaporated and the 
particles were sprayed over the probe volume. 

Different size particles could be introduced into the nebulizer thus 
producing monodisperse, btmodal, trimodal and quadrumodal distributions. The 
compressor was also used to trap the particles after they passed through the 
probe volume to avoid contaminating the surrounding environment. Figure 6 
shows a photograph of the optics and the particle generator. 

An alternative apparatus was also constructed to introduce glass beads 
into the probe volume. The glass beads come in dry containers and their size 
distribution was reasonably broad. A test tube with a stopper was used to 
contain the glass beads. Two glas3 tubes protruded through the stopper and 
flex tubing was connected at the end of each glass tube. One of the flex 
tubes was connected to a nitrogen tank via a pressure regulator and a flow 
meter. The other flex tube discharged the glass beads into the probe volume. 
The flow rate of nitrogen could be regulated very carefully and thus the flow 
velocity at the exit of the flexible tube could be reproduced. 


EXPERIMENTAL RESULTS 

The system was calibrated and tested with polystyrene latex particles of 
uniform and known size distributions. The size and uniformity of the 
particles were checked with a microscope and agreed with the manufacturer's 
specifications. 

Size and velocity distributions of the polystyrene particles flowing out 
of the heating chamber were obtained with the optical system described above. 
Figure 7a shows the distributions corresponding to 1.74 ym ai.u 3.3 yra. 

The calibration high voltage of 500V was established based the 1.74 ym. 

The arrow with the 3.3 ym mark indicates the value predicted by the Mie calcu- 
lations. In Figure 7b an intermediate size of 2.7 ym was added and the cali- 
bration high voltage was changed to 550V to decrease the size range. As 
before the calibration was based on the 1.74 ym, and the arrows point at the 
numerically predicted values. The 3.3 ym was measured very accurately, but 
the measurement of the 2.7 ym particles was off by two bins. This error was 
consistent and very repeatable. We attribute it to the oscillations found in 
the scattering function. Figure 7c shows the distributions corresponding to 
four different latex particles. As before the calibration was obtained with 
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Che 1.74 iim particles. The calibration high voltage is 550V. It tnust be 
pointed out that the size histograms are divided into 53 equal size bins. 
Therefore, distributions on the large diameter end will appear broader than 
those in the small diameter end. 

These results indicate that the technique and Instrument can accurately 
discriminate the particles by their size in LDV applications. Unfortunately, 
since all the particles are rather small and the flow velocity is also small 
and constant, all the particles moved at about the same velocity. Therefore, 
no velocity/size discrimination could be established. To illustrate this last 
point we used very large glass beads which would not follow the flow. Two 
different size classes were simultaneously introduced into the test tube and 
entrained into a constant nitrogen flow. These particles were sprayed over 
the probe volume and their size and velocity histograms are shown on Figure 8. 
Notice the bimodal size distribution and the corresponding blmodal velocity 
distribution. Although not obvious from this figure the small particles are 
moving faster than the big ones. This is better illustrated in Figure 9 which 
shows a velocity/size distribution of the entrained small and large glass 
beads. It can be deduced that the small glass beads are moving at about 
13 m/s while the large ones are moving at about 8 m/s. 

This particle size range is not expected in LDV applications but the 
experiment illustrates the kind of error that can be incurred in the flow 
velocity measurement if no consideration is given to the particle size. 

An aspect of the PIMAX technique that needs to be considered is the 
crosstalk between the two polarizations. To establish this crosstalk a series 
of tests were conducted. First, the polarization ratio of the transmitted 
laser beans was measured with and without particle extinction. To produce a 
substantial particle Interference a fuel spray was used to attenuate the laser 
beam by 10%. Very little change in the polarization ratio was measured which 
is of no surprise since the light scattered in the forward direction carries 
the same polarization as the incident radiation. A second set of tests was 
also conducted. Monodisperse droplets produced by a Berglund-Liu generator 
scattered light which could be analyzed by the receiving optics. One FMT 
measured the S-polarized light while the other measured the P-polarized light. 
We blocked the large beam which is S-polarlzed and illuminated the particles 
with the small P-polarized beams. The scattered light was measured in the PMT 
Intended for the S-polarized light. The crosstalk was about 4% at large (20°) 
angles and smaller at shallow angles. This crosstalk could be virtually 
eliminated by placing a polarizing filter in front of the PMT. It should be 
noted that the intensity into the small and large beams may not necessarily be 
the same and, therefore, the crosstalk will be proportional to the intensity 
ratio. For the small particles of interest in LDV and at shallow angles of 
collection this crosstalk is negligible. 

Obviously the two-color system (IMAX) is not subject to this crosstalk. 

In general, it is easier to separate the different colors of the scattered 
light than it is to separate the different polarizations. Yet budget con- 
straints may Impose the use of single color lasers in which case PIMAX is a 
viable technique. 
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CONCLUSIONS 

A method has been presented and results discussed to measure the size of 
particles i.» LDV applications. The accuracy and resolution of the method were 
illustrated by measuring polystyrene particles between 1.1 pm and 3.3 pm in 
diameter. The effect of crosstalk between the two polarizations was measured 
and established to be very small. 

It was shown, using large glass beads, that different size particles move 
at different velocities and, therefore, the velocity of the flow cannot be 
established from the particle velocity unless consideration to the particle 
size is given. 
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SIZE PARAMETER a 


Figure 2. Perpendicularly polarized scattered light 
intensity as a function of the size parameter a 
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(ft - 1.96 - 12.9 x 10 -7 ) 6 - 5°: ( ft - 1.66), 
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Figure 5. Schematic of PIMAX system 
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Figure 7c. Size and velocity histogram of 1.1 pm, 1.74 pm, 2.7 pm, and 3.3 

pm polystyrene spheres in air. Photomultiplier high voltage - 550 V. 
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Figure 9. v elocity distribution of small and large glass beads entrained 
in an air jet. 
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